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Abstract 



We propose a new cognitive framework for option price modelling, using quantum 
neural computation formalism. Briefly, when we apply a classical nonlinear neural- 
network learning to a linear quantum Schrodinger equation, as a result we get a nonlinear 
Schrodinger equation (NLS), performing as a quantum stochastic filter. In this paper, 
we present a bidirectional quantum associative memory model for the Black-Scholes-like 
option price evolution, consisting of a pair of coupled NLS equations, one governing the 
stochastic volatility and the other governing the option price, both self-organizing in an 
adaptive 'market heat potential', trained by continuous Hebbian learning. This stiff pair 
of NLS equations is numerically solved using the method of lines with adaptive step-size 
integrator. 

Keywords: Option price modelling. Quantum neural computation, nonlinear Schrodinger 
equations, leverage effect, bidirectional associative memory 

1 Introduction 

The celebrated Black-Scholes partial differential equation (PDE) describes the time-evolution 
of the market value of a stock option [H [2]. Formally, for a function u = u{t, s) defined on 
the domain 0<s<oo, 0<t<T and describing the market value of a stock option 
with the stock (asset) price s, the Black-Scholes PDE can be written (using the physicist 
notation: dzU = du/dz)as a diffusion-type equation0 which also resembles the backward 
Fokker-Planck equatioro 



Recall that this similarity with diffusion equation led Black and Scholes to obtain their option price formula 
as the solution of the diffusion equation with the initial and boundary conditions given by the option contract 
terms. 

^Recall that the forward Fokker-Planck equation (also known as the Kolmogorov forward equation, in which 
the probabilities diffuse outwards as time moves forwards) describes the time evolution of the probability 
density function p — p{t, x) for the position a; of a particle, and can be generalized to other observables as 
well [3]- Its first use was statistical description of Brownian motion of a particle in a fluid. Applied to the 
option-pricing process p — p{t,s) with drift Di = Di{t,s), diffusion D2 = D2{t,s) and volatility ct^, the 




(1) 



where cr > is the standard deviation, or volatility of s, r is the short-term prevailing 
continuously-compounded risk-free interest rate, and T > is the time to maturity of the 
stock option. In this formulation it is assumed that the underlying (typically the stock) 
follows a geometric Brownian motion with 'drift' fj, and volatility cr, given by the stochastic 
differential equation (SDE) [S] 

ds(t) = fj.s{t)dt + as{t)dW{t), (2) 

where W is the standard Wiener processH 

The solution of the PDE ([T|) depends on boundary conditions, subject to a number of in- 
terpretations, some requiring minor transformations of the basic BS equation or its solution]^ 
In practice, the volatility is the least known parameter in ([TJ, and its estimation is generally 
the most important part of pricing options. Usually, the volatility is given in a yearly basis, 
baselined to some standard, e.g., 252 trading days per year, or 360 or 365 calendar days. 
However, and especially after the 1987 crash, the geometric Brownian motion model and the 

forward Fokker-Planck equation reads: 

dtp = ^dss {D2(T^p) - ds {Dip) . 

The corresponding backward Fokker-Planck equation (which is probabilistic diffusion in reverse, i.e., starting 
at the final forecasts, the probabilities diffuse outwards as time moves backwards) reads: 

dtp = - ^cr^dss {D2p) - da (Dip) . 



^The economic ideas behind the Black-Scholes option pricing theory translated to the stochastic methods 
and concepts are as follows (see [7]). First, the option price depends on the stock price and this is a random 
variable evolving with time. Second, the efficient market hypothesis [8] [9], i.e., the market incorporates 
instantaneously any information concerning future market evolution, implies that the random term in the 
stochastic equation must be delta-correlated. That is: speculative prices are driven by white noise. It is 
known that any white noise can be written as a combination of the derivative of the Wiener process [34] and 
white shot noise (see [TD]). In this framework, the Black-Scholes option pricing method was first based on the 
geometric Brownian motion [Il[l], and it was lately extended to include white shot noise. 

The PDE IT} is usually derived from SDEs describing the geometric Brownian motion ([2} , with the solution 
given by: 

s(f) = s(0)e(''-^''''*+^^W. 

In mathematical finance, derivation is usually performed using Ito lemma [S^ (assuming that the underlying 
asset obeys the Ito SDE), while in econophysics it is performed using Stratonovich interpretation (assuming 
that the underlying asset obeys the Stratonovich SDE ^3]) |10l 17]. 

*The basic equation (|T]| can be applied to a number of one-dimensional models of interpretations of prices 
given to u, e.g., puts or calls, and to s, e.g., stocks or futures, dividends, etc. In the first (and most important) 
example, u{t, s) — c{t, s) is a call on a European vanilla option with exercise price X and maturity at T; then 
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BS formula were unable to reproduce the option price data of real markets O As an alter- 
native, models of financial dynamics based on two-dimensional diffusion processes, known as 
stochastic volatility (SV) models [11], are being widely accepted as a reasonable explanation 
for many empirical observations collected under the name of 'stylized facts' [12j. In such 
models the volatility, that is, the standard deviation of returns, originally thought to be a 
constant, is a random process coupled with the return in a SDE of the form similar to ([2]), 
so that they both form a two-dimensional diffusion process governed by a pair of Langevin 
equations [Tn[T3l[l4j. 

Using the standard Kolmogorov probability approach, instead of the market value of an 
option given by the Black-Scholes equation ([T|), we could consider the corresponding prob- 
ability density function (PDF) given by the backward Fokker-Planck equation (see |10j). 
Alternatively, we can obtain the same PDF (for the market value of a stock option), us- 
ing the quantum-probability formalism [39^ I40j . as a solution to a time-dependent linear 
Schrddinger equation for the evolution of the complex-valued wave ^/;— function for which the 
absolute square, IV'P, is the PDF (see [43J). 

In this paper we go a step further and propose a novel general quantum-probability 
basedi option-pricing model, which is both nonlinear (see [171 [HI [HI [20] ) and adaptive (see 

the solution to Q is given by (see, e.g. [6]) 

c(s,t) = siV(di) -Xe-''(^-''iV(d2), 
^ ln(g/X) + (r + la^)(r-t) 
' a(T-t)V2 

^ ln(./X) + (r-la^)(T-t) 
' a{T-ty/2 

where N{d) = i[l + Erf(d/A/2)]. 



^Recall that Black-Scholes model assumes that the underlying volatility is constant over the life of the 
derivative, and unaffected by the changes in the price level of the underlying. However, this model cannot 
explain long-observed features of the implied volatility surface such as volatility smile and skew, which indicate 
that implied volatility does tend to vary with respect to strike price and expiration. By assuming that the 
volatility of the underlying price is a stochastic process itself, rather than a constant, it becomes possible to 
model derivatives more accurately. 

®Note that the domain of validity of the 'quantum probability' is not restricted to the microscopic world |41) . 
There are macroscopic features of classically behaving systems, which cannot be explained without recourse 
to the quantum dynamics. For example, a field theoretic model leads to the view of the phase transition as 
a condensation that is comparable to the formation of fog and rain drops from water vapor. According to 
such a model, the production of activity with long-range correlation in complex systems like financial markets 
takes place through the mechanism of spontaneous breakdown of symmetry, which has been shown to describe 
long-range correlation in condensed matter physics. The adoption of such approach enables modelling of a 
financial markets and its hierarchy of components as a fully integrated macroscopic quantum system, namely 
as a macroscopic system which is a quantum system not in the trivial sense that it is made, like all existing 
matter, by quantum components such as atoms and molecules, but in the sense that some of its macroscopic 



3 



[21^ [6| I35j). More precisely, we propose a quantum neural computation [25j approach to 
option price modelling and simulation. Note that this approach is in spirit of our adaptive 
path integral [291 ED] approach to human cognition. 

2 The Model 

2.1 Bidirectional, spatio-temporal, complex-valued associative memory ma- 
chine 

The new model is defined as a self-organized system of two coupled nonlinear Schrodinger 
(NLS) equationsl3 one defining the option-price wave function tp = il;[t,s), with the corre- 
sponding option-price PDF defined by \ip{t,s)\'^, and the other defining the volatility wave 
function a = a{t,a), with the corresponding volatility PDF defined by \a{t,a)\'^. The two 
focusing NLS equations are coupled so that the volatility PDF is a parameter in the option- 
price NLS, while the option-price PDF is a parameter in the volatility NLS. In addition, both 
processes evolve in a common self-organizing market heat potential. 

Formally, we propose an adaptive, semi-symmetrically coupled, volatility -|- option-pricing 
model (with interest rate r, imaginary unit i = \/— T and Hebbian learning rate c), which 
represents a bidirectional NLS-based spatio-temporal associative memory. The model is 

properties can best be described with recourse to quantum dynamics (see [42] and references therein). 

^NLS can be viewed as a 'deformation' of the linear time-dependent Schrodinger equation from non- 
relativistic quantum mechanics. NLS no longer has a physical interpretation as the evolution of a quantum 
particle, but can be derived as a model for quantum media such as Bose-Einstein condensates (see e.g. \38\). 
In addition, for some values of parameters and initial/boundary conditions, the PDF-solution to NLS can 
be a soliton. This hypothetical option-price-soliton can potentially give a new light on financial markets 
simulation and will be explored elsewhere. 
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defined (in natural quantum units) by the following stiff NLS-systemU 

Volatility NLS : idta = -^s'^\'il;\^dssa + V{w)\a\'^a, (4) 

Option price NLS : idtip = — -s^jcrp + IV'I^V' + ^V'j (5) 

N 



with : V{w) = ''^^Wigi and 



i=l 



Adaptation ODE : Wi = —Wi + c\a\gi\ip\. (6) 

In the proposed model, the c-NLS (jl]) governs the short-range PDF-evolution for stochastic 
volatility, which plays the role of a nonlinear (variable) coefficient in ([5]); the -f/^-NLS ([5|) 
defines the long-range PDF-evolution for stock price, which plays the role of a nonlinear 
coefficient in dH. The purpose of this coupling is to generate a leverage effect, i.e. stock 
volatility is (negatively) correlated to stock return^ (see, e.g. |15j). The tu— ODE ([6]) defines 
the (cr, V')— based continuous Hebbian learning [361 |37]. The adaptive volatility potential 



*The proposed coupled NLS equations Q and ([5]) are both examples of the standard focusing NLS equation 
(see, e.g. gZ]): 

idt^ = -ldxx^ + V\i^f^, (3) 

where -i/) : R x R ^ C. It is a well-known fact that this equation has both soliton and multisoliton solutions (see 
[50ll49j ). Technically speaking, the NLS equation ((3)l is called 'mass and energy sub-critical' [48], 'scattering- 
critical' [ST], 'Galilean invariant' and 'completely integrable system' [521153] . arising from classical Hamiltonian 
energy function 



with the classical Poisson brackets 

{ip{x),ip{y)} = {ip*{x),ip*{y)} = 0, {ip*{x),ip{y)} = iS{s - y). 
Its quantum Hamiltonian operator has the form 



and the classical Poisson brackets are replaced by quantum Poisson brackets, or commutators: 

mx),i>{y)] = [r{x),i>'iy)] = 0, [i>'ix),i>{y)] = -Six - y). 
Linear component of ([3]l is the standard Schrodinger equation: 

idtTp = - ^dxx-ip + VTp. 

The NLS equation ([3]l has a 'blowup' whenever the Hamiltonian H is negative [541 [55] and in particular one 
has a blowup arbitrarily close to the ground state [561 157] . 

®The hypothesis that financial leverage can explain the leverage effect was first discussed by F. Black [31| . 
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V^w) is defined as a scalar product of the weiglit vector Wi and tlie Gaussian kernel vector 
Qi and can be related to the market temperature (which obeys Boltzman distribution [16j). 
The Gaussian vector gi is defined as: 

Qi = exp[-{d- niidf], d = ytarget-y, with (7) 
y = 2sin(60i), Vtarget = s\a\'^ds, {i = l,...,N), (8) 

where rrii = random(— 1.0, 1.0), while ds = {si — sq)/{N — 1) represents the stock-price 
increment defined using the method of lines (where each NLS was decomposed into N first- 
order ODEs; see Appendix). 

Note that each of the NLS equations ^ and ([5]) individually resembles 'quantum stochas- 
tic filtering/quantum neural network' models of |22[ 123^ 12^. Thus, the whole model effectively 
performs quantum neural computation [25j, by giving a spatio-temporal and quantum gener- 
alization of Kosko's BAM family of neural networks [26\ [27] . In addition, the solitary nature 
of NLS equations may describe brain-like effects frequently occurring in financial markets: 
volatility /price soliton propagation, reflection and collision (see [28]). 

2.2 Simulation results 

The model (jl])-(l6]) has been numerically solved for the following initial conditions (IC) and 
repeatable boundary conditions (BC): 

IC : cJi(0,s) = 0.25, V'i(0,s) = l, (9) 
BC : dta{t,so)=dta{t,si), d^t, sq) = O^t, si), (10) 

using the method of lines [44j (with = 30 lines per NLS discretization, see Appendix) 
within the Cash-Karp version [45] of the adaptive step-size Runge-Kutta-Fehlberg 4/5th 
order algorithm [36] . The average simulation time (depending on the random initial weights 
and Gaussians) was 10-30 seconds on a standard Pentium 4 PC, using Visual C+-|- compileJ^ 

A sample daily evolution of the volatility PDF is given in Figure [H The corresponding 
option-price PDF is given in Figure 121 with the Log-plot given in Figure |2j The corresponding 
Hebbian weights and Gaussian kernels are given in Figure HI 

2.3 Complex-valued wave-function lines 

Complex-valued quantum probability lines for option price evolution have some interesting 
properties of their own. Imaginary and real wave-functions corresponding to the option- 
price PDF from Figure [2] are given in Figure [H They can be combined in a phase-space 

^°In this respect it is important to note that in order to have a useful option pricing model, the speed of 
calculation is paramount. Especially for European options it is extremely important to have a fast calculation 
of the price and the hedge parameters, since these are the most liquidly traded financial options. 
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Figure 1: A sample daily evolution of the short-range volatility PDF: global surface (left) 
and individual lines (right). 

fashion to give complex-plane plots, as in Figure M These complex phase plots of individual 
wave-function lines depict small fluctuations around circular periodic motions. 

3 Conclusion 

In this paper we have proposed a new cognitive formalism for option price modelling, based on 
quantum neural computation. The model consists of a pair of coupled nonlinear Schrodinger 
equations, representing quantum-stochastic volatility and option-price evolutions. The two 
equations create a leverage effect by self-organizing in a common market-heat potential, 
using continuous Hebbian learning. An efficient numerical solver for this bidirectional spatio- 
temporal associative memory is implemented in C++ using the method of lines and adaptive 
step-size algorithm. Sample simulations are provided as well as the outline of the algorithm. 

4 Appendix: Method of Lines for Coupled NLS Eqs. (^-(^ 

Assuniing that we already have a stable adaptive Runge-Kutta-Fehlberg (RKF) integrator for systems 
of ODEs, we can implement the powerful method of lines [44] for the coupled NLS pair (jl])-® as 
follows. We start with the standard heat equation: 

dtip = ^d^xip, (11) 

where = ip{x,t) is a real-valued function. Solution of this PDF can be found if we have one initial 
condition: 

IC : ij{x,t = 0) = Vo (12) 
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Figure 2: A sample daily evolution of the long-range option-price PDF (corresponding to the 
volatility from Figure [T]): global surface (left) and individual lines (right). 



and two boundary conditions: 



BC : 'ip{x ~ xq, t) — ipij, where xq is the initial value of x, (13) 
dxtpi^ = x/, t) = 0, where x/ is the final value of x. 



A commonly used second order/central finite diff'erence approximation for is 



Aa;2 

where k is an index designating a position along a grid in x which has M (e.g. 10) points and Ax is 
the spacing in x along the grid. In this way, the heat PDE (jlip can be approximated as a system of 
ODEs 

^ (fc ^ 1. 2, . . . TV = number of ODEs). (14) 

dt 2 Ax^ ^ ' ^ ' 

The set of ODEs ([H)) is then integrated, using the RKF integrator, subject to IC and BC. 

Once we are able accurately solve the set of ODEs with the conditions (|12p - P^ . we can add 
the potential term V{x): 

dti) = -dxxi' + V{x) tp, 
with the corresponding set of approximating ODEs: 

^4^^^±^^^|i^ + nfc)^., (. = 1,2,... TV). (15) 

Once we are able to solve the set of heat ODEs (fT5|) with the real- valued conditions (fT2| ~ (|13p . 
we can move into the complex plane, by introducing the imaginary unit i = \/— T on the left and the 
minus sign on the right: 

idtiP = -l-dx.4' + V{x)iP, (16) 




Figure 3: Log-plot of the long-range option-price PDF from Figure [2j 

This is the linear time-dependent Schrodinger equation from non-relativistic quantum mechanics (in 
natural units, for which both the Planck constant and the particle mass are one, see |30 ) for the 
complex-valued wave ^/;— function. It has the corresponding set of approximating complex-valued 
ODEs: 

'^^-2 A^^ + V{k)^,, {k = l,2,...N). (17) 

Once we are able to solve the linear Schrodinger equation (|16l) with the complex-valued conditions 
similar to P^ - p^ . as a set of complex- valued ODEs ([T7]), we can replace the linear term on the 
right with the cubic nonlinearity: 

idt^j = -^d.,^tP + V{xM^^. (18) 

This is the nonlinear Schrodinger equation (NLS). Its corresponding set of approximating complex- 
valued ODEs reads: 



^ _± v..+i , ^ vix.M'i^, (fc = 1, 2, . . . N). (19) 



, dipk ^ 1 V-'fc+i - 2V'fc + i^k-i 
dt 2 Aa;2 

Finally, our option pricing model (m)-® represents a stiff pair of coupled complex-valued ODE- 
decompositions of the form of (fT9| . together with the Hebbian learning ODE ([6]). Although the 
resulting ODE-system is stiff, the Cash-Karp version [551 of the adaptive step-size Runge-Kutta- 
Fehlberg algorithm [46] can efficiently solve it (in Visual C++) for various boundary conditions and 
random initial weights and Gaussian kernels. 
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